
TECH LIBRARY KAFB, NM 





TECH LIBRARY KAFB, NM 


■iiinii 

01M30n 


REPORT 1280 


THEORETICAL INVESTIGATION OF FLUTTER OF 
TWO-DIMENSIONAL FLAT PANELS WITH ONE 
SURFACE EXPOSED TO SUPERSONIC 
POTENTIAL FLOW 

By HERBERT C. NELSON and HERBERT J. CUNNINGHAM 


Langley Aeronautical Laboratory 
Langley Field, Va. 


National Advisory Committee for Aeronautics 

Headquarters, 1512 H Street NW., Washington 25, D. G. 

Created by act of Congress approved March 3, 1915, for the supervision and direction of the scientific study 
of the problems of flight (U. S. Code, title 50, sec. 151). Its membership was increased from 12 to 15 by act 
approved March 2, 1929, and to 17 by act approved May 25, 1948. The members are appointed by the President, 
and serve as such without compensation. 


Jerome C. Hunsaker, Sc. D., Massachusetts Institute of Technology, Chairman 
Leonard Carmichael, Ph. D., Secretary, Smithsonian Institution, Vice Chairman 


Joseph P. Adams, LL.B., Vice Chairman, Civil Aeronautics 
Board. 

Allen V. Astin, Ph. D., Director, National Bureau of Standards. 
Preston R. Bassett, M. A., Vice President, Sperry Rand Corp. 
Detlev W. Bronk, Ph. IX, President, Rockefeller Institute for 
Medical Research. 

Thomas S. Combs, Vice Admiral, United States Navy, Deputy 
Chief of Naval Operations (Air). 

Frederick C. Crawford, So. IX, Chairman of the Board, 
Thompson Products, Inc. 

James H. Doolittle, Sc. D., Vice President, Shell Oil Co. 
Clifford C. Furnas, Ph. D., Assistant Secretary of Defense 
(Research and Development) Department of Defense 


Carl J. Pfingstag, Rear Admiral, United States Navy, Assistant 
Chief for Field Activities, Bureau of Aeronautics. 

Donald L. Putt, Lieutenant General, United States Air Force, 
Deputy Chief of Staff (Development). 

Arthur E. Raymond, Sc. D., Vice President — Engineering, 
Douglas Aircraft Co., Inc. 

Francis W. Reichelderfer, Sc. D„ Chief, United States 
Weather Bureau. 

Edward V. Rickenbacker, Sc. D., Chairman of the Board, 
Eastern Air Lines, Inc. 

Louis S. Rothschild, Ph. B., Under Secretary of Commerce 
for Transportation. 

Nathan F. Twining, General, United States Air Force, Chief 
of Staff. 


Hugh L. Dryden, Ph. D., Director 


John F. Victory, LL. D., Executive Secretary 


John W. Crowley, Jr., B. S., Associate Director for Research Edward H. Chamberlin, Executive Officer 


Henry J. E. Reid, D. Eng., Director, Langley Aeronautical Laboratory, Langley Field, Va. 
Smith J. DeFrance, D. Eng., Director, Ames Aeronautical Laboratory, Moffett Field, Calif. 
Edward R. Sharp, Sc. D., Director, Lewis Flight Propulsion Laboratory, Cleveland, Ohio 

Walter C. Williams, B. S., Chief, High-Speed Flight Station, Edwards, Calif. 


n 


REPORT 1280 


THEORETICAL INVESTIGATION OF FLUTTER OF TWO-DIMENSIONAL FLAT PANELS WITH 
ONE SURFACE EXPOSED TO SUPERSONIC POTENTIAL FLOW 1 

By Herbert C. Nelson and Herbert J. Cunningham 


SUMMARY 

A Rayleigh type analysis involving chosen modes of the 
panel as degrees of freedom is used to treat the flutter of a two- 
dimensional flat panel supported, at its leading and trailing 
edges and subjected, to a middle-plane tensile force. The 
panel has a supersonic stream passing over its upper surface 
and still air below. The aerodynamic forces due to the super- 
sonic stream, are obtained from the theory for linearised 
two-dimensional unsteady flow and the forces due to the still 
air are obtained from acoustical theory. 

In order to study the effect, of increasing the number of modes 
in the analysis , two and then jour modes are employed. The 
modes used are the first four natural modes of the pane! in a 
vacuum with no tensile force acting. The analysis includes 
these variables: Mach number, structural damping, tensile 

force, density of the still air, and edge fixity ( clamped and 
pinned). For certain combinations of these, variables, stability 
boundaries are obtained which, can be used to determine the 
panel thickness required to prevent flutter for any panel material 
and altitude. 

In contrast to some previous panel-flutter investigations, 
the results of the present analysis show that, sufficiently thick 
panels are flutter free for the Mach numbers treated and suggest 
that this is true throughout the supersonic speed range. 

A comparison of results from the present theory for flat 
panels and from a criterion developed by II. P. Isaacs for the 
static stability of buckled panels is made with a few experimental 
results on flat and buckled jmnels clamped a,t leading and 
trailing edges. 

INTRODUCTION 

The flutter of thin metal plates or panels, such as com- 
pose the covei'ing or skin of missiles and other craft intended 
for high-speed flight, has recently become of increased con- 
cern. Such panels may be initially flat or curved and may 
be small or fairly large in aspect ratio. In addition, they 
may be prestressed and will probably become warped in 
flight by aerodynamic heating. If one or more of the 
panels on a particular configuration are vibrating, tiie 
basic structure supporting them can usually be considered 
rigid. The fixity at the edges of the panels ranges between 
clamped and pinned, depending on the construction. Some 
preliminary experimentation and analytical work suggests 
that this type of instability is of concern only at supersonic 
speeds. 


The problem of panel flutter embraces so many possible 
factors as to discourage general treatment, and previous 
papers on the subject (for example, refs. 1 to 7) have em- 
ployed various simplifying assumptions in order to obtain 
specific solutions to what might perhaps be considered 
different, phases of the problem. In all the references cited, 
(lie main assumption made is that a panel and the flow over 
it are two-dimensional. Other assumptions common to 
the reference papers are that small-deflection plate theory 
and linearized How theory may be used. 

References 1 lo 4 examine tin' ease of a panel buckled by a 
constant shortening and held at its leading and trailing edges, 
with a supersonic stream over its upper surface and no per- 
turbation pressures on its lower surface. In reference 1 
steady-state air forces and in reference 2 quasi-stationary air 
forces, which include the first order of the frequency of oscil- 
lation, are used. Both these references consider the dynamic 
stability of the buckled panel. Reference 3 and the more ex- 
haustive reference 4 examine the static stability of t lie 
buckled panel and propose that motion (flutter) is t lie result 
when static equilibrium is not possible. Reference 5 and a 
section of reference 2 treat the ease of a flat panel pinned at 
its leading and trailing edges. Reference 5 uses exact linear- 
ized unsteady aerodynamic forces and therefore, in contrast 
to reference 2, imposes no limitations on the order of the fre- 
quency. In references 1, 2, and 5 a generalized-coordinate 
approach involving chosen modes of the panel as degrees of 
freedom is employed. Reference 6, on the other hand, indi- 
cates how the problem of a vibrating membrane in a super- 
sonic stream can be treated bv means of Laplace transforms 
and suggests that similar treatment can be given to the plate 
problem. Reference 7 treats the ease of a two-dimensional 
panel on many equally spaced simple supports with compres- 
sible air flowing over the upper surface and dead air below 
the panel, and the results indicate that the possible panel 
instabilities are divergence for subsonic flow and flutter for 
supersonic flow. Some questionable features of the results 
obtained in references 2, 5, and 7 are examined in the section 
entitled “Results and Discussion” in the present report. 

A Rayleigh type flutter analysis is developed herein by 
means of Galerkin’s process for a two-dimensional flat panel 
held in some manner at its leading and trailing edges and 
acted on by a middle-plane or axial force (which, in the case 
of tension, introduces a restoring force similar to that for the 
membrane). The upper surface of the panel is subjected to 


1 Supersedes NACA Technical Note 34G5 by Herbert C. Nelson and Herbert J. Cunningham, 1955. 
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a supersonic stream and the lower surface to an unconfined 
mass of stationary compressible air. The normal modes of 
the panel with no middle-plane force acting are used as de- 
grees of freedom in the anaNsis. As in reference 5, exact 
linearized unsteady aerodynamic forces are employed. In 
the reference paper the integrals yielding these forces are 
evaluated analytically. In the present report these integrals 
are evaluated numerically. 

Numerical results are presented in order to examine some 
effects of including two and then four modes in the analysis 
and to determine effects of Mach number, density of the 
supersonic stream, panel mass and stiffness, edge fixity (to 
some extent), structural damping, axial load, and density of 
the still air below the panel. In appendix A an alternative 
solution by means of Laplace transforms is developed for the 
plate problem just described. No numerical results are ob- 
tained by this met hod, however. 

A comparison of results from the theory presented herein 
for flat panels and from a criterion of reference 3 for the static 
stability of buckled panels is made with a few experimental 
results for flat and buckled panels clamped at tin* leading 
and trailing edges. 
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SYMBOLS 

speed of sound in undisturbed medium 
structural and aerodynamic integrals de- 
fined after equation (13) 
panel chord 

Et 3 

local flexural stiffness, , .r N 

’ 12(1 — v 1 ) 

Young’s modulus of elasticity of panel 
material 

tension parameter, F/c’ 2 m A o> i 2 
functions defined in equation (24) 
external force per unit width acting in 
midplane of panel (tensile force posi- 
tive) 

structural damping coefficient 
matrix elements defined in equation (13) 
Hank el function of second kind, of zero 
order, (notation of ref. 18) 
aerodynamic integrals defined after equa- 
tion (2.3) 

Bessel functions of order p of first and 
second kind, respectively, (notation of 
ref. 18) 

reduced frequency, cu/2l ' 
stiffness parameter (reduced first natural 
frequency), cccJ2l’ 

eigenvalues defin'd after equations (16) 
and (17) and given for first four panel 
modes in table I 

aerodynamic functions defined after equa- 
tion (21) 

aerodynamic functions defined after equa- 
tion (Blfi) 


m A panel mass per unit surface area, <rr 

M Mach number, U/a 

N u No coefficients in mode-shape equations (16) 

and (17), given for first four panel 
modes in table I 

p(x,t) net perturbation pressure, positive down- 

ward 

p„( x) pressure coefficient associated with mode 

shape Z„, defined in equation (20) 
p,„Pi upper- and lower-surface contributions to 

perturbation pressure, respectively 
]>,„,)>„ pressures in undisturbed supersonic stream 

_ and still-air region, respectively 

y > „(.r),/ , „(.r),/k,('') components of p„{x) defined after equa- 
tion (22) 

<1 dynamic pressure, ^ pi' 1 
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coefficients of equations (28) and (Bio) 
tabulated in appendix B for first, four 
modes of panel with clamped edges 
time 

velocity of supersonic stream 
coordinates defined in figure 1 
vertical displacement of panel 
flutter mode shape 

wth natural mode shape for panel vibrating 
in vacuo 

panel-air mass ratio, m A /pc 
Poisson’s ratio 

densities in undisturbed supersonic stream 
and still-air region, respectively 
density of panel material 
local thickness of panel 
disturbance-velocity potential 
frequency of oscillation 
frequency associated with mode shape Z„ 
frequency parameter, 2kM 2 /p 2 

frequency ratio squared \ , except in 

flutter calculations where P = ) (1 Fig) 

square matrix- 
column matrix 


Subscripts: 

u upper-surface contribution 

I lower-surface contribution 

Primes denote differentiation with respect to the argument. 


ANALYSIS 

STATEMENT OF PROBLEM 

A thin isotropic, two-dimensional plate (beam) of con- 
stant thickness, as shown in figure 1, is considered herein. 
The plate, is undergoing simple harmonic motion and is 
acted on by a middle-plane or axial force V (tension or 
compression); its upper surface is subjected to a supersonic 
stream of velocity U, pressure and density p, and its 
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lower surface is subjected to still air with pressure p„ and 
density p 0 . The differential equation of motion for tin- 
plate may be written as 

D dr ‘ +m " 1 

where, in the present case, the vertical displacement of tin- 
plate z(r,t) may be expressed as Z(x)e tut , w is the circular 
frequency of oscillation, is the net perturbation pres- 

sure (positive downward) arising from the motion of tin- 
plate, m A is the plate mass per unit surface area, and the 
local flexural stiffness D is given by /£r 3 /12( 1— r 2 ). For the 
case where p(x,t) = 0, equation (1) may be obtained, with 
appropriate changes in notation, from reference 8. 

In the remaining development the constant-pressure 
term — of equation (1) is considered to be zero. This 
in no way affects the generality of the results for the oscil- 
lating plate, since inclusion of the constant-pressure term 
as nonzero would result only in adding a particular solution 
which represents a. static vertical deflection. In addition, 
the coordinate ,r of equation (1) is divided by the plate 
chord c and henceforth is employed in this nondimensional 
sense. Thus, equation (1) multiplied by e~ ul becomes 

D t Z"" — trm A Z~Z"-\-p(.r J 1)e-‘ u, ==0 (2) 


appendix A are the solutions to these problems by means 
of Laplace transforms. 

NET PERTURBATION PRESSURE p(x,t) 

The net pressure p{x,t), as mentioned previously, arises 
from the oscillatory motion of the plate. It is this pressure 
which damps, or in the case of flutter sustains, the oscillation. 
The pressure on the upper surface is obtained from the 
theory for linearized unsteady supersonic flow and the 
pressure on the lower surface from acoustical theory. The 
perturbation pressure in terms of the pressures p u on the 
upper surface and pi on the lower surface is 



p(x,t)=p„—pi 

(5) 

where 

/ d<j> u [I d<A\ 

(0) 


p, ‘~~ p \dt + c bx) 

and 

b4>t 

Pl = pa d t 

(7) 


From reference 9 the velocity potential for the upper surface 
can be obtained in the form 



Based on reference 10, the velocity potential for the lower 
surface can be obtained, as shown in appendix B, in the form 




where //„ <2> (.r) is the Hankcl function of (he second kind, 
of zero order. 


SOLUTION BY GALERKIN’S METHOD 


where the primes denote differentiation with respect to 
the argument x. 

In order to obtain a specific solution of equation (2), 
four boundary conditions are required. The plate is con- 
sidered to be held at its leading and trailing edges as shown 
in figure 1, and this assumption leads to the conditions for 
pinned edges: 

Z(0)=Z(1) = Z"(0)=Z"(1) = 0 (3) 

and for clamped edges: 

Z(())-Z(l)=Z , (0)-Z'(l)-0 (4) 

In a later section of the analysis the boundary-value prob- 
lems, as exemplified by equations (2) and (3) or (2) and 
(4), are solved by Galerkin’s method. Also considered in 


Outline of method. — The boundary-value problems con- 
sidered earlier (eqs. (2) and (3) for the pinned-edge plate 
and equations (2) and (4) for the clamped-edge plate) are 
now solved by means of Galerkin’s method. (A detailed 
account of Galcrkin’s method may be found in ref. 11.) As 
a first step, the flutter mode shape Z(x) is approximated by 
a linear combination of the form 

Z(x)=cj>^,a„Z„(x) (10) 

n - 1 

where the coefficients a„ may represent complex amplitudes 
and where the functions Z„(.r) are the mode shapes for the 
plate vibrating in a vacuum without an axial force /'’acting. 
The function Z t is the fundamental mode shape associated 
with the lowest natural frequency coj , and the remaining 
functions Z.., Z : „ . . . Z N are consecutively the higher modes 
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shapes. The shapes Z n satisfy the pinned-edge or clamped- 
edgc boundary conditions {Z n replacing Z in eqs. (3) and (4)) 
and the differential equation 

1 ^Z,:'"-^m A Z„=() ( 11 ) 


Thus the flutter condition may be expressed in the form 


Gn a* • • • Gw 
G n G», ... G m 


=0 


( 14 ) 


where is the frequency of oscillation for which Z„ is the 
mode shape. 

The remainder of the tlalerkin process for solving the 
aforementioned boundary-value problems consists of deter- 
mining the coefficients a„ of equation (10) in the following 
manner: Substitute equation (10) into equation (2), replace 

the term ^Z,,"" by <*i„‘m A Z„ in accordance with equation 

(11), multiply by one of the mode shapes Z„, integrate the 
result from r=0 to x—1, and equate to zero. When n is 
made 1 , 2, . . . N in succession, N linear equations are 
obtained which determine the unknowns a„. These equa- 
tions can be written in the form 


G n 

G» ■ 

■ 0i.v 


'a, ' 


'O' 



Go i 

Goo . 

• G-2N 

-< 

a-2 

► = H 
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The matrix elements are given by 

G m „ = ll-^A m „ — Q j J 1 ('mil (13) 

where 


m. A 

M== 

pc 

A m „= 

fz,„Z„ (lx 
/ 0 



r l z,„'z n ’dx 

In 

f~ = ~2 /Y 2 

c J m A 

== 

' m ft 

( Z m p„(x)d.r 
J 0 


and where p„{x) is the. pressure )>(x,t), obtained from equa- 
tions (5) to (9) with Z replaced by Z„, multiplied by 
e~ twl lpcco 2 . In equation (13) a), and a>„ are the first and the 
mth natural frequency, respectively, of the plate with no 
axial force acting. 

Flutter determinant.— Tlie flutter condition or condition 
of harmonic vibration, which is given by the nontrivial 
solution for the coefficients a„, is obtained from equation (12) 
by setting the determinant of the matrix G,„ „ equal to zero. 


Gn I G N 2 ... (7; 


■*NN 


Remarks on alternative procedure. — The procedure from 
equation (11) to equation (13) is, in general, not the most 
accurate that could he followed for values of F other than 
zero. A generally more accurate procedure would be to use, 
instead of equation (11), the differential equation for the 
panel with tension: 


^ Z n ""-^m A Z„-^Z n "= 0 


(15) 


When equation (15) is solved, subject to the appropriate 
boundary conditions, the frequencies co„ are found to be 
functions of F for both pinned and clamped edges, but the 
mode shapes Z, t do not vary with F for pinned edges. The 
use of equation (15) rather than equation (11) would mean 
that in equation (13) the term JB„ ul would not appear and 
the frequencies and mode shapes would be those that satisfy 
equation (15). 

Equation (11) rather than equation (15) has been used 
herein for the following reasons: For pinned-edge panels 
there is no difference in the mode shapes or in the final 
numerical flutter results; for clamped-edge panels the deter- 
mination of the values of Z„ and w„ that satisfy equation 
(15) is laborious and must be carried out for every desired 
value of F. Elimination of the term //)„,„ from the matrix 
elements, through use of equation (15), does not compensate 
for the labor of determining the natural frequencies and mode 
shapes as functions of F. The differences in final numerical 
flutter results for the clamped-edge panel approach zero as 
the number of modes in the analysis is increased and are 
expected to be small even when only a few modes are used. 

EVALUATION OF TERMS IN FLUTTER DETERMINANT 


Structural integrals A mn and B mn and frequencies w n . — 
Consideration will now be given to the evaluation of the 
mode-shape integrals and frequencies in the elements of 
equation (15). The mode shapes Z„ and associated natural 
frequencies u n obtained from equation (12) are: 

For the pinned-edge plate, 


Z„ = N, sin K„.r 

u„—K„ 2 I P 

\ fn A r ^ 


( 10 ) 


THEORETICAL INVESTIGATION OF FLUTTER OF FLAT PANELS EXPOSED TO SUPERSONIC POTENTIAL FLOW 


5 


where K n is obtained from the frequency equation 


sin K „= 0 


For tho clamped-edge plate, 

Z„=A ? 1 [cos K n x — cosh jK^z-f-A^Csin K n x — si] ih K„x)] 


D 

m A c 4 


y (i7) 


TABLE I— MODE-SHAPE FACTORS, EIGENVALUES, 
FREQUENCY RATIOS, AND STRUCTURAL INTEGRALS 
FOR FIRST FOUR NORMAL MODES 

(a) Pinncd-cdge plate * 


Mode, n 

JV, 

N 2 

K n 

Wn/ui 

Ann 

Bnn 

1 

1 


X 

1 

0.5 

0. 5x- 

2 

1 


2x 

4 

.5 

2. Ox 2 

3 

1 


3ir 

0 

.5 

4. Sir 2 

4 

1 


4 x 

16 

.5 

8. Ox 2 


*A m n=Bmn = Q 


where K„ is obtained from the frequency equation 
cos K„ cosh K„ = 1 

In equations (16) and (17) the factor ATi is used to produce 
unit deflection at the center of the plate (z=0.5) for modes 
that are symmetric about the center and at the point of 
maximum deflection between the leading edge (x=0) and 
the center of tho plate for modes that are antisymmetric 
about the center. The factor N-> in equation (17) is estab- 
lished by the boundary condition requiring zero deflection 
at, x= 1 and is expressed by 

(18) 

sin K n — sink K n 

The quantities w It co„,/c a,, A mn , and B mn required in equa- 
tion (14) can be determined directly from equations (16) or 
(17). First, however, values must be established for N 2 and 
K n in the case of the pinned-edge plate and for N u N 2 , and 
K„ in the case of the clamped-edge plate. Table I includes 
values of all these quantities for the first four modes of 
vibration. 

The values for A mn shown in table I arc zero when m^n 
because of the orthogonality of the mode shapes Z„ of equa- 
tions (16) and (17). For the pinned-edge case the slopes 


(b) Clamped-edge plate b * 0 


Mode, 

n 

Ni 

Ni 

K n 

lOnfal 1 

A„n 

Bnn 

1 

2 

3 

4 

-0. 629690 
-.66260 
. 7109645 
-. 06120074 

-0. 98250 
-1.000778 
-.99997 
-1.00000145 

4.730 

7. 853204625 
10. 99560784 
14. 13716549 

1 

2. 7566 
5. 404 
8. 933 

0. 396 
.440 
.506 
.432 

4. 88 
21.2 
49.8 
76.4 


A mn =0 (w/n) 

« i3 =B 3 i = 4.302 

7*24= -8.300 

Other limn—Q (niF^u) 

« The significant figures shown were used to avoid sin all -difference errors in mode shapes 

Z,' of the mode shapes arc also orthogonal and, consequently, 
B mn is zero when mZn. For the clamped-edge case, oven 
though the slopes Z„' are not orthogonal, the integrand of 
B mn is antisymmetric about x=0.5 when m and n are not 
both even or both odd and, consequently, B nn is zero when 
m^n except for B vi , B su B 2i , and B l2 (for the first four 
modes). 

Aerodynamic integrals Cmn. The remaining term in the 

elements of equation (14) that requires evaluation is the 
integral 

Cm"=j o Z m p„{x)d.r, (19) 

As mentioned previously, p„(x) is the pressure p(x,t), obtained 
from equations (5) to (9) with Z replaced by Z, n multiplied 
by e~ iat /pca> 2 . The quantity p„(x) is therefore given by 


ir [- 2 -®+i f ] [h <*-»] '« } ).+ 

7 [2 £ ZMHo m (2kM\x-it\)d^ ( 20 ) 

where the contributions from the upper and lower surfaces of the plate are designated by subscripts u and l, respectively. 
Upon elimination of the derivatives in the integrands of the upper-surface contribution through integration by parts, 
performance of the indicated differentiation of the second integral, and extraction of the singularity at £=x in the 
Hankel function of the lower-surface contribution, equation (20) may be written in the form 


Pn(x) = 


Iff 

/3 LJo 


Z„(mx-m+-~ M p- Z„(x)+^ Z/(r)] M 


+ P5 


[j: 


if 

7T Jo 


Z n (£)L(x — £)(/£+- Z„(£) logc (kM\x- 


Tl)^] ; 


( 21 ) 
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where 



, .2 M 

+* T' Jl 



L\v)= 


' J</2kM\u\)+[ 


Y n (2kM\u\)- 


.log, (kM\u\) 


A(0)= 7 +: 7=0.57721 (> (Killer’s constant) 

7 r 1 


The quantities J„ and Y „ are t he Bessel functions of order 
p, of the first and second kind, respectively. 

For convenience p„{ jc) is considered in three parts, namely, 

P„ (-r) - l\ (r ) + T\ (.r) + F„ (.r) (22) 


where 


pm=\ [JJ 2 z/(r)] 

p„u)=- f 1 z n {S)U*-m 

P Jo 

p„U)- Po 1 [ ' z,m log, (kM\.r-a\)da 

P IT Jo 


Hence, equation (19) may be put in the form 


where 


(dm ti I m a I I m >i I J mv 


/„,„ = J o Z„l\{x)<Ir 

/,„„=( Z,„P„{.r)d.r 


7 




■/. 


ZJ\,(x)ds 


(23) 


Tlie first integral represents the effect of the supersonic 
stream passing over the upper surface of the plate; the other 
two integrals represent the effect of the still air below the 
plate. 

Before further development of the method of the present 
report for determining the aerodynamic treatments of 

references 1, 2, and 5, which deal with the pinned-edge 
plate, will be examined. These references consider only the 
effects of the supersonic stream, the air below the plate 
being treated, in essence, as massless; that is, p 0 /p is taken to 


be zero and the integrals I mn and I mn are omitted. In 
references 1 and 2 the aerodynamic effects are accounted for 
as if the integral I mn has been expanded as a power series in 
the frequency of oscillation; reference 1 retains only the 
steady-state or zero-order frequency term and reference 2 
adds the first-order frequency term. In reference 5, on the 
other hand, the integral is evaluated exactly with regard 
to the frequency. This is possible because the modal 
functions Z„ for the pinned-edge plate are sine waves (see 
eq. (1C)) so that can be obtained in terms of the functions 
(sometimes called Schwarz functions) 

fJ«,f»)= x 1 M ( ( U ) (hi (X = 0 , 1 ) (24) 

h J I) \<t / 

where 



O) 


and b has the four values 

*-=+{!:}' 

A similar result could be obtained for the damped-edge 
plate by approximating the modal functions Z n (see eq. (17)) 
by a finite sine series 

Ft 

Z „ « T'.tL sin rir.r. 

r = 1 

For either pinned or clamped edges, the arguments a and b 
of the Schwarz functions f\ would range from large positive 
to large negative values, particularly for the higher modes, 
and would thus require extensive tabulation of f Q and A. 
The exact expression for the pressure term />„(■>') is employed 
in the present report but. because the necessary tables of / 0 and 
/, are not available, for convenience, a numerical method of 
integration is used to evaluate p„(x) and the aerodynamic, 
integrals and of equation (23). 

The numerical method is based on the following integration 
rules for parabolic arcs: 

£ ’ y(x)dx=^ (5j/(ri)+8i/(r 2 )— y(r 3 )] (25a) 

j J y(j')dx=Y^ [— yi-ri)Jr8y(-r 2 )-\-5y(x 3 )] (25b) 

where = .ri -|- Ar and / 3 =/ 2 +A/. The range of integration 
in equation (19), 0^/^l, is, for convenience, divided into 
an even number of equal segments. Fi-om the standpoint 
of accuracy the number of segments needed depends on the 
number of nodes in the highest mode and on the value of 
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to for which p„(x) is evaluated. For the mimerical applications of the present report, 10 segments were found to be 
adequate, and the method of integration is illustrated for this number of segments in the equations to follow. 

The use of 10 segments would, in general, require the determination of p n (x) in equation (19) at 11 points on the plate. 
However, the integrand of C m „ at x=0 and x=l is zero since the mode shapes Z m are zero at these points, and therefore 
p n (x ) need be evaluated at only the 9 interior points (equally spaced between x=0 and x=l). The values of the terms 
P n , P n , and P„ of equation (22) for these points may be arranged in matrix form as follows : 


>»(.!)' 


“ 82/(0) 

— 2/( — . 

P«(. 2) 


132/(.l) 

72/(0) 

Pn(.Z) 


132(.2) 

122/(.l) 

Pn(A) 

1 

*" 1 20/3 

132/(.3) 

122/(.2) 

- Pn(. 5) 

132/(.4) 

122/(.3) 

Pn(- 6 ) 

132/(.5) 

122/(.4) 

Pn(- 7) 


132/ (.6) 

122/(.5) 

Pn(. 8) 


132/(.7) 

122/(.6) 

_ Pn(.9) , 


_132(.8) 

122/(.7) 


0 

0 

0 

0 

— 2 (— . 

1) 0 

0 

0 

72,(0) 

-2-(-.I) 

0 

0 

122(.l) 

72(0) 

— L{—. 

1) 0 

122/(.2) 

122(.l) 

72(0) 

-2(— .1) 

12L (.3) 

112(.2) 

152(.l) 

42(0) 

122/(.4) 

112(.3) 

142(.2) 

122(.l) 

122(.5) 

112(.4) 

142(,3) 

112(,2) 

122(.6) 

1 12(.5) 

142(.4) 

1 12(.3) 


0 

0 

o” 


r ^„(.i) 

0 

0 

0 


X„(.2) 

0 

0 

0 


Z„(.3) 

0 

0 

0 


'/■n(A) 

0 

0 

0 


Z,.(.5) 

0 

0 

0 


X„(.6) 

52(0) 

0 

0 


*n(.7) 

132(.l) 

52(0) 

0 


X„(.8) 

122(.2) 

132(. 1) 

52(0)_ 


. Z„(.9) 





x„(.2) 


Z,/(.2) 

Z n (.3) 


Z„'(.S) 

Z„(. 4) 


Z„'(A) 

" Z„(-5) 

Z n '(.5) 

Z„(.6) 

Z„'l 6) 

*.(.7) 



X„(.8) 


X„'(.8) 

^n(.9) > 


^/(•9)J 


>»(.!)' 


2(0) 

2(.l) 

2(.2) 

2(.3) 

2(.4) 

2(.5) 

2(.6) 

2(.7) 

2(.8) 


' 13^„(.l)' 

Pn{.2) 


2(.l) 

2(0) 

2(.l) 

2 (.2) 

2(.3) 

2(.4) 

2(.5) 

2(.6) 

2(.7) 


\27. n (.2) 

Pn(. 3) 


2 (.2) 

2(.l) 

2(0) 

2(. 1) 

2(.2) 

2(.3) 

2(.4) 

2(. 5) 

2(.6) 


12Z n L3) 

P n(-4) 

Po/p 

120 

2(.3) 

2(.2) 

2(.l) 

2(0) 

2(. 1) 

2(.2) 

2(.3) 

2(.4) 

2(.5) 


11Z„(A) 

- 2*„(.5) 

2 (.4) 

2(.3) 

2(.2) 

2(.l) 

2(0) 

2(.l) 

2(.2) 

2(.3) 

2(.4) 

-< 

14Z„(.5) >- 

2„(.6) 

2(.5) 

2(.4) 

2 (.3) 

2(.2) 

2(. 1) 

2(0) 

2(.l) 

2(.2) 

2(.3) 


1 1X„(.6) 

2„(.7) 


2(.6) 

2 (.5) 

2(.4) 

2(.3) 

2(.2) 

2(. 1) 

2(0) 

2(. I) 

2(.2) 


I2X„(.7) 

Pn(-S) 


2(.7) 

2(.6) 

2(.5) 

2(.4) 

2 (.3) 

2(.2) 

2(.l) 

2(0) 

2(.l) 


12Z„(.8) 

„ P /i(-9) _ 


_2(.8) 

2 (.7) 

2(-6) 

2(.5) 

2(.4) 

2(,3) 

2(.2) 

2(. 1) 

2(0) _ 


. 13Z„(.9) _ 


’>»(.l)' 1 


2i(.l) 

2o(.l) 

2 3 (.l) 

2<(.l) 

2s(.l) 

2«(.l) 

2-(.l) 

2 8 (.l) 

2„(.l) 

2, 0 (.l) 


Si 

P„(.2) 


2, (.2) 

2,(.2) 

2 3 (.2) 

2d. 2) 

2 S (.2) 

2 e (.2) 

2, (.2) 

2 8 (.2) 

2„(.2) 

2io(.2) 


s 2 

P„(.3) 


2i(.3) 

2, (.3) 

2,i (.3) 

2d-3) 

2s(.3) 

Zs(.3) 

Lr(.3) 

2„(.3) 

2 9 (.3) 

2, 0 (.3) 


s 3 

Pn{A) 


2, (.4) 

2 2 (.4) 

7/3 (-4) 

La (.4) 

2 5 (.4) 

2e(.4) 

Li(A) 

2s(.4) 

2,(.4) 

2, 0 (.4) 


S 4 

„ Pn(- 5) 

, po/p 

2i(.5) 

2 2 (.5) 

2 3 (.5) 

2d-5) 

2 S (.5) 

L 6 (.5) 

Lit 5) 

2 8 (.5) 

2s(.5) 

2,„(.5) 

-< 

Ss - 

Fn(-6) 

4 

2i(.6) 

2 2 (.6) 

2 3 (. 6) 

Zc (.6) 

2 5 (.6) 

Let 6) 

Iv(.6) 

2,(.6) 

2„(.6) 

2io(-6) 


Se 

P n(*7) 


2, (.7) 

2 2 (.7) 

2s(.7) 

2c (.7) 

2s (-7) 

2a(.7) 

2t(.7) 

2s(.7) 

2o(.7) 

2jo(.7) 


s? 

Pn(-8) 


2, (.8) 

2 2 (.8) 

2 3 (.8) 

2c(-8) 

2 5 (.8) 

2»(.8) 

Li ( * 8) 

2 g (.8) 

2 S (.8) 

2, 0 (.8) 


s 8 

P„(. 9) 


2i(.9) 

2 2 (.9) 

2 3 (.9) 

2j(-9) 

2s(.9) 

2«(.9) 

L 7 (.9? 

2,(.9) 

2„(.9) 

2io(-9) 


So 












_ 




(26) 


(27) 


(28) 


where L(u) and L(u) are defined after equation (21) and L r (u) and s T are defined in appendix B (eqs. (B16) and (B15)). 
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The rows of the square matrix in equation (26), down to the 
row pertaining to P n (. 5), were obtained by applying the 
integration rule given by equation (25a). In the remaining 
rows the contributions from the region z>0.5 were obtained 
by applying equation (25b). In equation (27) the inte- 
grating factors multiplying Z n in the column matrix were 
obtained by using equation (25a) between x—Q and £=0.5 
and equation (25b) between z=0.5 and x=1.0. In appendix 

B the singular integral P n (%) presented after equation (22) 

is evaluated and the form of P n leading to equation (28) is 
derived. 

Equations (26), (27), and (28) arc summed in accordance 
with equation (22) to obtain the column matrix {p„(xj\. 
By use of this column matrix and the integration rules of 
equation (25), the aerodynamic term C mn is obtained in 
the form 

<7 m „= J ^ [1 3p„(.l)Z m (.l)+12p„(.2)Z OT (.2)+12p B (.3)Z m (.3)+ 

llp„(.4)Z m (.4)+14p„(.5)Z M (.5) + llp„(.6)Z m (.6) + 

1 2p n (. 7 )Z m ( . 7) + 1 2p„(. 8)Z„,( . 8) + 1 3p„(. 9)Z m ( . 9)] (29) 

where the integrating factors 13, 12, . . . 12, 13 were ob- 
tained in the same manner as those in equation (27). By 
means of equation (29), C mn can be evaluated for a given edge 
fixity and for particular values of M, k, and p 0 /p. 

SOLUTION OF FLUTTER DETERMINANT 

As previously stated, the conditions for flutter are de- 
termined from the nontrivial solutions of equation (14). 
Since equation (14) is complex, it may be solved directly for 
one complex unknown or two real unknowns. For a specific 
edge fixity the variables (see eq. (13)) in equation (14) are 
1/p (the inverse of p is preferred because p becomes infinite 
for p=0), if, /, M, k, and po/p. It is convenient to interpret 
the 9 of equation (13) as the complex quantity (&j,/cj) 2 (1 +ig) 
rather than (W“) 2 , where g may be regarded as a structural 
damping coefficient. (For this use of g, sec, for example, 
refs. 12 and 13.) Each of the various quantities on which 
equation (14) is dependent was varied to some extent, as 
will be discussed in the next section. A particular calcula- 
tion was performed by setting values for 1/p, M, k, /, and 
p 0 /p and solving for 9. Then, because it was one of the more 
easily varied parameters, 1/p was changed and again Q was 
solved for. This procedure was continued until curves of 
1/p and ( R.P.ay * 2 plotted against g passed through g = 0. 
The value for k was then changed and the procedure re- 
peated. After sufficient variation of 1/p and k, curves 
could be established of 1/p against 2k l = 2k{R.P.9l) U2 for 
particular values of the other parameters M, g, f, and p 0 /p. 

RESULTS AND DISCUSSION 

In the preceding sections a method of flutter analysis has 
been developed for a two-dimensional flat panel or plate held 
at its leading and trailing edges. The variables in the 


analysis are the number of modes or degrees of freedom the 
panel is assumed to have, Mach number (greater than 1.0), 
1/p, 2k 1 =2k(R.P.Q) 112 , g, f, po/p, and edge fixity. The 
analysis conveniently yields stability boundaries in terms of 
1/p and 2k i , which are used as the coordinates of most of the 
figures presented. These two parameters are given in terms 
of the properties of the panel and supersonic stream by 

1_ p c 

1 ^ KS r j p E [ (30) 

" 1 U V24(l -V 2 ) c \cg, 

where cr is panel density, q is dynamic pressure, and K, is the 
first-mode eigenvalue given in table I for clamped and pinned 
edges. Inasmuch as the various parameters in the analysis 
contain implicitly the panel properties (E, <r, v, and r/c), 
axial force F, air density, and speed of sound, the effects of 
varying these implicit properties can be obtained only by 
cross-plotting. 

Some effects of the number of modes used in the analysis 
are studied by using two and then four modes of the clamped- 
edge panel with selected values of M, g, /, and p 0 /p. In 
addition, M, g, /, and p 0 /p ai'e varied in order to study their 
effects. To a lesser extent the pinned-edge panel is investi- 
gated for comparison with certain clamped-edge results. 

The following table lists the conditions for which stability 
boundaries are given : 


Mach 

number, 

M 

Degrees 

of 

freedom 

Structural damping 
coefficient, g 

Tension param- 
eter, / 

Density 

ratio, 

PqJp 

Clamped edge 

1 

0 

0 

0 

1.3 


0, 0.1, 0.5, 1.0 

1 . 0 

4 

0, 0.005, 0.025, 0.03, 0.05 

0 

0 

0 

0, 0.1, 0.5, 1.0 

0 

72 

2 

0, 0.002, 0.00375, 0.05 

0 

0 

1.56 

2 

0 

0 

0 

4 

0 

0 

0 

Pinned edge 

72 

2 

0. 0.003, 0.00475, 0.05 0 

0 


The results are first grouped according to Mach number 
and are later summarized and compared. 

RESULTS FOR MACH NUMBER OF 1.3 

Effects of two and four modes. — Figure 2 gives the results 
for the clamped-edge panel for the simple case of two degrees 
of freedom (first and second modes) with g~f—p 0 /p=0. 
The abscissa is the stiffness parameter co i c/U=2k l and the 
ordinate is the mass ratio 1/p. An ordinate of zero repre- 
sents the limiting case of p=(), or, in other words, the plate 
is vibrating in a vacuum. The two solid curves are the first- 
and second-mode stability boundaries as indicated. It was 
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Stiffness parameter, U| c/U 


Figure 2. — Stability boundaries from a two-mode analysis for 
clamped-edge panels. M= 1.3; g=f=—- =0. 

p 

established by application of the Nyquist criterion (see, for 
example, ref. 14) as well as by interpretation of structural- 
damping results that the region to the right of the first mode 
boundary is stable, whereas the region to the left is unstable; 
furthermore, the region within the second-mode boundary is 
doubly unstable as indicated (unstable with regard to both 
boundaries). Values of the reduced frequency k are indi- 
cated along both curves. The points at which the curves 
cross the abscissa correspond to vibration in a pure normal 
mode (flutter at the limiting condition of p= 0). 

It can be seen from equations (30), by taking the product 
of 1 In and 2 k u that a specified panel material, air density, 
and speed of sound are represented by a hyperbola such as 
the dashed curve of figure 2 with the panel thickness-chord 
ratio r/c increasing to the right. The intersection of the 
hyperbola with the stability boundary fixes the value of r/c 
for neutral stability. Thicker panels are stable and thinner 
panels are unstable. (The particular hyperbola shown is for 
aluminum panels in air with standard sea-level properties. 
For denser panels or less dense air, the hyperbola would be 
below the one shown.) 

Some effects of the number of modes in the analysis were 


studied by including the first four normal modes, and the 
results are shown as the solid curves of figure 3. The dashed 
curves are the results for two modes from figure 2. With the 
addition of the third and fourth modes, the first-mode bound- 
ary is moved very slightly to the left (except where it crosses 
the abscissa) and is still the “critical” or decisive stability 
boundary separating the stable from the unstable region. 
The second-mode boundary is also only moderately affected. 
Within the already unstable region there now exist third- 
mode and fourth-mode boundaries which are closely anal- 
ogous in appearance and character to the first- and second- 
mode boundaries, respectively. The unstable region is 
divided by three of the boundaries, into regions of different 
degrees of instability as indicated by the numbers in paren- 
theses ranging from (1) to (4). (The points at which the 



Figure 3. — Stability boundaries from two-mode and four-mode 
analyses for clamped-edge panels. M=1.3; g=f=~= 0. 
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various boundaries cross the abscissa have the same signifi- 
cance as before.) These results indicate that two modes 
give a decisive stability boundary which is a close approxima- 
tion to that for a large number of modes, at least for the con- 
ditions g=f—polp=0 and M= 1.3. 

Effects of structural damping coefficient g.- — Figures 4 (a) 
and 4 (b) show the first-mode and second-mode boundaries 
(from a four-mode analysis) for various values of g (taken 
to be the same for all modes). Third- and fourth-mode 
boundaries are affected by g in a manner similar to that of 
the first and second modes, respectively, and are not shown. 
The second-mode boundary of figure 4 (b) vanishes com- 
pletely when g becomes slightly greater than 0.025, and for 
all positive values of g it remains in the unstable region to the 
left of the first-mode boundary. Included in figure 4 is the 
dashed hyperbola from figure 2. Since the thickness-chord 
ratio t/c decreases to the left in the figure, the abscissas of 
the intersections of the hyperbola with the stability bound- 
aries in figure 4 (a) show the proportional reduction in 
thickness required to prevent flutter as g increases. 

Effects of tension. — Tension has a marked effect on the 



O .2 .4 .6 .8 1.0 


Stiffness parameter, W| c/U 
(a) First-mode boundary. 

Figure 4. — First-mode and second-mode stability boundaries from 
a four-mode analysis for clamped-edge panels for various values of 
structural damping coefficient g. M= 1.3; /=^=0. 


stability boundaries, as shown in figures 5 (a) and 5 (b). 
Figure 5 (a) shows the pertinent segments of the first-, 
second-, third-, and fourth-mode stability boundaries for 
g = 0 and for the three values 0.1, 0.5, and 1.0 of the tension 
parameter/. As/ increases, all the boundaries move to the 
left, and the thickness required to prevent flutter is de- 
creased. Furthermore, as/ increases, the first-mode bound- 
ary moves to the left more rapidly than the higher mode 
boundaries so that the rightmost boundary, separating 
stable from unstable regions, is one of the higher mode 
boundaries. For example, for /= 1.0 in figure 5 (a), the 
third-mode boundary is farthest to the right. This trend is 
not surprising since application of tension to the clamped- 
edge plate causes the largest percentage increase in the first 
natural frequency, the next largest in the second natural 
frequency, and so on. Thus, it appears that the inclusion 
of only two modes in a flutter an alysis may not be sufficient 
when the plate is subjected to tension. Inasmuch as the 
stiffness parameter 2 k, and the tension parameter / arc both 



O .1 .2 .3 .4 .5 .6 .7 

Stiffness parameter, u, c/U 


(b) Second-mode boundary. 
Figure 4. — Concluded. 
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0 .1 .2 .3 .4 ,5 .6 .7 

Stiffness parameter, 


(a) 5 = 0. 

(b) g = 0.005. 

Figure 5. — Stability boundaries from a four-mode analysis for 
clamped-edge panels for three values of the tension parameter /. 

M =1.3; ^=0. 


based on the first natural frequency of the panel without 
tension, the shift of the boundaries is due solely to the 
tensile force F. 

Figure 5 (b) shows segments of the first-, second-, third-, 
and fourth-mode boundaries with j — 0. 1 , 0.5, and 1.0 for 
g=0.005. By comparing figures 5 (b) and 5 (a) it can be 
seen that g has a marked effect for the smaller values of / but 
its effect diminishes as,/ increases. 

Effects of still air below panel. — The one remaining pa- 
rameter to be considered at M= 1.3 is p 0 /p, the ratio of the 
density of the still air below the panel to the density of the 
supersonic stream above. Tn the preceding results this 
ratio was zero. The effect of increasing p 0 / p to 1.0 will now 
be examined. For the sake of simplicity and convenience, 
only a two-mode analysis is made. Effects of structural 
damping and tension are also included. 

Figure 6 (a) shows first- and second-mode boundaries for 
p 0 /p= 1 as solid curves and, for comparison, the dashed 
boundaries for p 0 /p= 0 from figure 2. Just as with the other 
results, the points where the boundaries cross the abscissa 


correspond to pure-mode resonance in a vacuum. At these 
crossings the imaginary part of C„„ passes through zero. 
This imaginary part is a measure of aerodynamic damping. 
In the previous calculations C nn consisted only of 
whereas, for p 0 /p=1.0, C nn also contains I„ n -\-T nn (see eq. 
(23)). By comparison of the solid and dashed curves on 
figure 6 (a) it can be seen that, as a consequence, the first- 
mode boundary has moved to the left by about 20 percent 
but the second-mode boundary has changed relatively little. 

Such an effect of still air might be expected since, for the 
same maximum panel amplitude, a first-mode vibration 
radiates into the still-air region a greater amount of energy 
per cycle than does a second-mode vibration. (With regard 
to the radiation of sound from a piston in a plane wall, 
specifically for the case of a piston with nonrigid face, p. 336 
of ref. 15 gives the result that, at frequencies which are 
small compared with the ratio of the speed of sound to 2ir 
times the piston radius, the pressure on the piston is approxi- 
mately uniform and nearly proportional to the average 
velocity of the piston. Since the average velocity of the 
second mode and all other antisymmetric modes is zero, the 
pressure due to these modes is nearly zero and, accordingly, 
almost no work is being done on the still air.) From the 
fact that net energy can never pass from the still air into the 
panel, it does not follow, however, that the still air neces- 
sarily has a stabilizing effect in all cases. Conceivably, the 
still air could act to modify the flutter mode so that more 
energy would be extracted from the supersonic stream, and 
thus contribute toward an instability. Apparently such is 
the case in figure 6 (a), where the solid second-mode curve 
is above the dashed second-mode curve. The fact that 
dissipation of energy into the still air is not necessarily 
stabilizing should not be surprising, inasmuch as another 
means of energy dissipation, structural damping, is usually 
stabilizing but sometimes destabilizing. 

As can be observed in figure 6 (a), the first-mode boundary 
has moved to the left of the second-mode boundary in the 
region of small mass ratio; in this region the second-mode 
boundary becomes critical. 

Figure 6 (b) shows the effects of structural damping on the 
first-mode boundary, which for g—0 is shown more com- 
pletely in figure 6 (a). Curves are included for g=0, 0.005, 
0.03, and 0.05. For values of g larger than about 0.025 the 
second-mode boundary vanishes as it did previously with 
Po /p= 0 in figure 4 (b), and only the first-mode boundary 
remains. The dashed hyperbola for aluminum at sea level 
is included in figure 6 (b) , and it can be seen that a plate with 
zero structural damping would have to be about 30 percent 
thicker than one with g=0.05 in order to prevent flutter. 

Figure 7 shows the effects of the tension parameter / for 
Po /p= i and g = 0. Both first- and second-mode boundaries 
are shown for /=0, 0.1, 0.5, and 1.0. In this case, just as 
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(a) First-mode and second-mode boundaries for 0. 

(b) First-mode boundary for various values of g. 

Figure 6. — Stability boundaries from a two-mode analysis for 
clamped-edge panels for various values of g. M= 1.3; /=0; — = 0 
and 1.0. 

with po/p— 0, tension causes a marked reduction in the thick- 
ness required to prevent flutter. Furthermore, if more than 
two modes had been included, tension would have resulted 
in a higher mode boundary farther to the right than the 
curves shown for the higher values of /. 

RESULTS FOR MACH NUMBER OF ^ 

Clamped-edge panels. — Figure 8 (a) shows the stability 
boundaries obtained from a two-mode analysis for clamped- 
edge panels at with/=p 0 /p=0 for various values of g. 

Included in the figure is the dashed hyperbola appropriate 
to this Mach number for aluminum panels in sea-level air. 
From a comparison of figure 2 and the curves of figure 8 (a) 



Figure 7. — Stability boundaries from a two-mode analysis for 
clamped-edge panels for various values of/. lW=1.3:g = 0; —=1.0 

p 


for g=0, the first-mode boundary of figure 2 appears to have 
moved into the positive mass-ratio region and the second- 
mode boundary appears to be moving toward the negative 
mass-ratio region. Such is the case, but, inasmuch as the 
flutter frequencies on the upper boundary of figure 8 (a) are 
about midway between the first and second natural fre- 
quencies, this boundary can now be referred to only loosely 
as a “first-mode” boundary. The lower boundary is still 
readily identified as a second-mode boundary and the inter- 
section with the abscissa corresponds to vibration in a pure 
second mode. 

In contrast to the situation at M= 1.3, the second-mode 
boundary for ^ = 0 is now decisive for panels represented by 
the dashed hyperbola. Values of thickness to the right of the 
second-mode boundary are stable and, in addition, a small 
range of thickness values is stable between the upper and 
second-mode boundaries. 

The curves in figure 8 (a) for positive values of g show that 
the region of instability within the second-mode boundary 
is reduced for small values of g (as for M— 1.3) and vanishes 
when g is slightly greater than 0.00375, but that small values 
of g increase the region of instability associated with the 
upper boundary. This effect of g on the upper boundary is 
in marked contrast to its effect on the first-mode boundary 
at M= 1.3. (See fig. 4 (a).) The differing effects of struc- 
tural damping at M= 1.3, M=- v / 2, and M=1.56 are con- 
sidered further in the section on “Variations With Mach 
Number.” 

Pinned-edge panels.— In order to indicate effects of edge 
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fixity, boundaries are shown in figure 8 (b) for conditions 
identical to those of figure 8 (a) except that the edges are 
pinned rather than clamped. The boundaries for y — 0, 
which are given incompletely in reference 5 and thereby lead 
to the conclusion that only a small range of panel thickness 
is stable at have been extended to higher frequencies 

with the result that sufficiently thick panels are also found 
to be stable. The effect of structural damping on both 
boundaries in figure 8 (b) is very similar to that in figure 8 (a) . 
The dashed hyperbola appropriate to pinned-edge aluminum 
panels in sea-level air is included in the figure. The hyper- 
bolas of figures 8 (a) and 8 (b) are located differently because 
of the different values for the first-mode eigenvalue K x of 
equation (30) for pinned and clamped edges. (See table I.) 
From the intersection of the dashed hyperbolas with the 
stability boundaries in the two figures, it can be determined 
that a pinned-edge panel must be somewhat thicker than a 
clamped-eclge panel in order to be flutter free but not nearly 
as thick as might be expected from a simple comparison of 
the first natural frequencies. Values of the reduced fre- 



0 .1 .2 .3 4 .5 .6 .7 .8 


Stiffness parameter, a^c/U 

(a) Clamped-edge panel. 

(b) Pinned-edge panel. 

Figure 8. — Stability boundaries from a two-mode analysis for 
clamped-edge and pinned-edge panels for various values of g. 

jtf=V2;/=^°=0. 

p 


quency k are indicated along each of the boundaries of 
figure 8. 

Based on what occurred at M= 1.3 (see fig. 3), there is the 
possibility that for g = 0 the fourth-mode boundary from a 
four-mode analysis would alter the stability picture in both 
figures 8 (a) and 8 (b) in the relatively unimportant narrow 
range of stability between the boundaries shown for g=0. 
This minor effect of the fourth-mode boundary is expected 
to disappear for values of g greater than about 0.005 and, 
therefore, a four-mode analysis was not made for this Mach 
number. 

RESULTS FOR MACH NUMBER OF 1.56 

Effects of two and four modes. — As in tbe case of M= 1.3, 
stability boundaries were obtained first for two and then for 
four degrees of freedom with g=/=po/p=0. These bound- 
aries appear in figure 9 as dashed curves for two modes and 
solid curves for four modes. Values of the reduced frequency 
k are indicated along the boundaries. The stable region is 
again to the right, and on the left the degree of instability 
is indicated in parentheses for the four-mode analysis. 

The two-mode results in figure 9 continue the trend noted 
in the preceding section from comparison of the curves of 
figure 2 and those of figure 8 (a) for£ = 0. The second-mode 



Figure 9. — Stability boundaries from two-mode and four-mode 

analyses for clamped-edge panels. M= 1.56; g=f=— = 0. 

P 



14 


REPORT 1280 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


boundary has moved entirely into the negative mass-ratio 
region. The upper boundary has moved higher in the posi- 
tive mass-ratio region, and the flutter frequencies along it, 
which in figure 8 (a) were midway between the first and 
second natural frequencies, are now closer to the second. 
For this reason the upper boundary, which was loosely identi- 
fied as a “first-mode” boundary in the discussion of figure 
8 (a) , will now be referred to as a “second-mode” boundary. 

A further point of difference between the results at M— 1 .3 
and M=1.56 is that the addition of the third and fourth 
modes at M= 1.56 shifts the decisive stability boundary to 
the left by about 10 percent, whereas at M=1.3 the shift is 
insignificant. (Compare figs. 3 and 9.) Although this shift 
indicates that the two-mode result is not well converged, the 
two-mode boundary is conservative ; that is, it requires a 
greater thickness to prevent flutter. (As with two modes, 
when four modes are used, half of the stability boundaries fall 
in the negative mass-ratio region.) 

Effects of structural damping coefficient g. — No curves are 
shown to indicate effects of structural damping at a Mach 
number of 1.56, the reason being that, for moderate values 
of the coefficient g, ranging at least up to 0.05, the stability 
boundaries fall virtually on top of those for^ = 0. The major 
effect of structural damping is a moderate change in flutter 
frequency. 

Effects of tension and of still air below panel. — Effects of 
tension have not been determined, but tension is expected 
to have essentially the same favorable stiffening effect at all 
Mach numbers as at M— 1.3. The effect of still air behind 
the panel lias also not been determined, but this effect is 
expected to be less than at M= 1.3 for two reasons: First, the 
air beneath the panel acts primarily as an energy absorber 
and one means of energy absorption, structural damping, 
has been found ineffective in shifting the stability boundaries. 
Second, on the decisive boundary the flutter mode appears to 
be predominantly the second natural panel mode, and it was 
found that at M= 1.3 the second-mode boundary is changed 
only slightly by increasing p 0 /p from 0 to 1. 

VARIATIONS WITH MACH NUMBER 

The foregoing results have been presented for particular 
Mach numbers. In an effort to clarify some of the anom- 
alies that have been noted in these results, figures 10 to 12 
are presented. Figure 10 shows the panel thickness-chord 
ratio required to prevent flutter as a function of M for 
clamped-edge panels with g=f=p 0 /p=0. The curves apply 
to aluminum panels in standard sea-level air. The values 
at M= 1.3, y'2, and 1.56 were obtained from figures 2 and 
8 (a) and the two-mode results of figure 9. The shape of the 
curves between these known points is estimated. The 
stable region is above or to the right of the shaded boundaries. 

The boundaiy which is labeled “first-mode” on one end 
and “second-mode” on the other has flutter frequencies 
which progress from slightly above the first natural frequency 
to somewhat below the second natural frequency as the 
Mach number is increased. (See previous discussions con- 
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Figure 10. — Minimum panel thickness ratio r/c required to prevent 
flutter as a function of Mach number for clamped-edge aluminum 
panels in sea-level air. g=f=— =0. 


cerning figs. 2, 8 (a) , and 9.) The boundary labeled “second- 
mode” has flutter frequencies slightly below the second 
natural frequency throughout. 

Figure 10 shows the second-mode stability boundary to 
be decisive in the Mach number range from slightly above 
1.30 to slightly above -y/2. As the structural damping g 
is increased from zero, the second-mode boundary shrinks 
to the left leaving the “first-mode” — “second-mode” bound- 
ary decisive throughout the range of M shown. For example, 
for a value of g slightly greater than 0.0038 the second-mode 
boundary does not exist at M=-yj 2 (see fig. 8 (a)), and for 
a value of g slightly greater than 0.025 it does not exist at 
M— 1.3 (see fig. 4 (b)). 

These effects of g on the second-mode boundary are 
illustrated in figure 11, which contains cross plots of g 
against r/c obtained from the intersections of the dashed 
hyperbolas (for aluminum panels in sea-level air) with the 
boundaries for constant g such as shown in figures 4, 8 (a), 
and 9. Figure 11 also shows that an increase in g from 
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(a) M = 1.30. 

(b) M=V 2. 

(c) M= 1.56. 

Figure 11. — Structural damping coefficient g as a function of panel 
thickness ratio t/c from a two-mode analysis for clamped-edge 

aluminum panels in sea-level air. f= — = 0. (Refer to dashed 

p 

hyperbolas in figs. 2, 8 (a), and 0.) 



Figure 12.— Minimum panel thickness ratio r/c required to prevent 
flutter as a function of Mach number for aluminum panels in sea- 
level air. 


zero would cause the “first-mode” — “second-mode” bound- 
ary of figure 10 to drop markedly at M= 1.3, rise slightly 
at M=- y'2, and remain essentially unchanged at M= 1.56. 
The ratios of flutter frequency to the first natural frequency 
co/coi are indicated for each of the crossings and tend to 
show more clearly the connection between figures 10 and 11. 
(Values of w/m near w 2 /oh, which is approximately 2.76, are 
associated with the second-mode boundary and values 
between 1.0 and 2.05 are associated with the “first-mode” — 
“second-mode” boundary.) A complete understanding of 
the manner in which the curves change character and position 
with Mach number, particularly between M— 1.3 and 
requires more calculation than presented herein. 

Figure 12, which has the same ordinates as figure 10, is 
presented for the purpose of summarizing some effects of 
all the parameters investigated. The results shown are 
based on two modes, except in the case of tension where only 
four-mode results are known. The results again apply to 
aluminum panels in sea-level air. The figure shows as a 
solid curve the shaded boundary from figure 10 for clamped- 
edge panels and as a short-dash curve the effect on this 
boundary of increasing g from 0 to 0.05. The third (long- 
dash) curve is for pinned-edge panels with y = 0, the value 
at M= 2 having been obtained from reference 2. The points 
at M=-\/ 2 were obtained from figure 8 (b) and the upper 
(second-mode) curve was patterned after that for clamped- 
edge panels. As a matter of interest, points are included 
in figure 12 at M=1.3 for clamped-edge panels and indicate 
the effects of tension (/= 0.5) and of still air below the panel 
(po/p=1.0) forg=0 and £ = 0.05. 

Some effects of the various parameters can be assessed 
from figure 12. The overall result is that r/c is highest in 
the low supersonic Mach number range and suggests that 
this range is the more critical from a design standpoint. 
Structural damping is seen to have a large favorable effect 
near and below M=y' 2. Although rather influential at 
M= 1.3, the still air below the panel is expected to have less 
effect at A/=V2 and 1.56. Tension, which is seen to have a 
large favorable effect at M— 1.3, is expected to be similarly 
effective for all Mach numbers. In this connection, it might 
be mentioned that one means of producing tension is by a 
static-pressure difference between the upper and lower panel 
surfaces, particularly for the case where the panel leading 
and trailing edges are prevented from moving toward each 
other. A comparison of the results for the edge fixities, 
pinned and clamped, is of interest because the edge fixity of 
actual panels falls somewhere between. 

COMPARISON WITH OTHER THEORETICAL WORK 

In reference 2 the conclusion is reached that all panels, 
regardless of thickness, are unstable for supersonic Mach 
numbers less than ^2. This result and the more plausible 
results of reference 2 for are based on air forces 

expanded to the first power of the frequency of oscillation. 
In reference 5 the necessity of including higher order fre- 
quency terms for Mach numbers near y/2 is pointed out, and 
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stability boundaries, based on exact linearized unsteady air 
forces, are presented for M—-J2 and M= 2. One bound- 
ary is obtained at M= 2 which agrees well with the com- 
parable result from reference 2, whereas two boundaries are 
obtained at M=-*j2. The boundaries at M=- y/2, because 
they are not carried high enough in frequency, are interpreted 
in reference 5 as showing that stability is possible only for a 
small range of panel thickness at this Mach number. The 
results of reference 2 for M<C\/ 2 are not questioned in 
reference 5. 

In the present report, stability boundaries are computed 
for M=^j2 and for Mach numbers above and below this 
value (namely, M= 1.3 and M— 1.56). In view of the find- 
ings of references 2 and 5, perhaps the most noteworthy 
result of the present investigation is that, for the Mach 
numbers treated and probably throughout the supersonic 
range, sufficiently thick panels ax - e stable. 

In references 2 and 5 and the present report, 
appears as a transitional value. The transition is evidenced 
herein by the contrasting behavior of the stability boundaries 
at M= 1.3 and M= 1.56. Some understanding of why a 
Mach number of ^/2 is transitional can be had by examining- 

matrix equation (26). The term ^ {Z„} of equation 


(26), being the entire first-order frequency contribution to 
the damping, is dominant at low frequencies. This term 
appears to control the slope, at low frequencies, of the 
eventually decisive stability boundary and changes sign as M 
passes through A /2. When M <C y'2, the slope is negative for 
low frequencies, but as the frequency increases the slope 
eventually becomes positive because of the higher order 
frequency effects (for example, in fig. 2). Because only 
first-order frequency effects are included, in essence only the 
beginning portions of the stability boundaries for MO-^/2 are 
obtained in reference 2, and, as a consequence, the conclusion 
is reached that all panels are unstable below this Mach 
number. For M> A /2 the slope of the decisive stability 
boundary starts out positive and becomes more so as the 
frequency increases. (See fig. 9.) If aspect ratio were 
included in the present treatment (by considering three- 
dimensional rather than two-dimensional panels), a reduction 
in aspect ratio would probably tend to eliminate the initial 
negative slope of the eventually decisive stability boundary 
for MO V2 and increase the initial positive slope for 
This effect of aspect ratio is expected because, in general, a 
reduction in aspect ratio results in an increase in aerodynamic 
damping with a consequent enlargement of regions of 
stability. 

In reference 7, which treats a different problem (namely, an 
infinite two-dimensional panel on equally spaced supports), 
the result was also obtained that somewhere in the supersonic 
Mach number range a panel will flutter regardless of its 
thickness. The conclusion was reached that stability is not 
possible at supersonic Mach numbers less than about 1.25 
and that at higher Mach numbers a sufficient increase in 
thickness will always render a stable panel unstable. How- 
ever, it was observed that over a large portion of the pre- 
dicted region of instability the flutter was of an extremely 


mild character, since a large number of oscillations were 
required to double the amplitude. With the hope of elimi- 
nating the large region of mild instability, small amounts of 
viscous damping were included. Contrary to expectations, 
thick panels remained unstable for the example given at 
M=1.8. 

As part of the viscous-damping investigation, the results 
were interpreted so as to determine regions of stability and 
instability. As shown in figure 10 of reference 7, an apparent 
conflict with the results of Nyquist diagrams was found. 
(The Nyquist diagram concept is used in general in reference 
7 for determining stability.) This conflict is based on the 
assumption that most investigators interpret structural- 
damping results according to the concept that removal of 
damping tends to destabilize. This assumption is incorrect, 
however, and no such simple criterion holds true for inter- 
preting structural-damping results. A feature to be noted 
in the example chosen in reference 7 to illustrate the apparent 
conflict is the existence of infinite singularities in the air forces 
at the end points of the boundaries (g=0, 0.01, and 0.03) on 
the right in figure 9 of the reference (designated type B loci 
therein). By way of explanation, such singularities occur hi 
the Iinearized-fiow treatment because a traveling wave of 
panel deflection is moving at a speed corresponding to M= 1 
relative to the air above or below the panel. 

In this analysis the question of stability was investigated 
by means of both the structural-damping concept and the 
Nyquist diagram concept. The structural-damping results 
in every case agreed with the Nyquist diagram results. 
Incidentally, in using the Nyquist concept, knowledge of the 
aerodynamic forces for all frequencies from minus infinity 
to plus infinity is required. Thus, the concept is not appli- 
cable, in general, when the air forces are approximated by a 
few terms of a power-series expansion in the frequency of 
oscillation. 

COMPARISON WITH EXPERIMENT 

A few experimental results on the flutter of flat and buckled 
panels are available for comparison with the theory of the 
present report for flat panels and that of reference 3 for 
buckled panels. Reference 16 gives experimental results 
at M= 1 .3 for panels 11.62 niches long in the stream direction 
and 8 inches wide that were held by clamping the leading and 
trailing edges. In figure 13, the results of reference 16 at 
M— 1.3, together with data more recently obtained in the 
Langley supersonic flutter apparatus on both flat and 
buckled panels at other Mach numbers, are compared with 
theory. The results arc presented in terms of the thickness- 
chord ratio t/c needed to prevent flutter of aluminum-alloy 
panels at an altitude of 25,000 feet as a function of Mach 
number. These points were obtained from tests of panels 
of different thicknesses (see, for example, ref. 16) and repre- 
sent the thinnest panels which did not flutter. (Where 
necessaiy, experimental data were adjusted to a pressure al- 
titude of 25,000 feet with the relation t/c— (t/c) r (q/q r ) U3 . 
The subscript r refers to the experimental conditions.) 

In figure 13, the solid curve is the flutter boundary for 
flat panels obtained from the present theory and the square 
symbols are the corresponding experimental results. The 
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Figure 13. — Minimum panel thickness ratio r/c required to prevent 
flutter of aluminum-alloy panels at 25,000-ft altitude. 

clashed curve is the static-stability boundary for buckled 
panels, proposed as a flutter boundary in reference 3, and the 
circles are the associated experimental results. The theo- 
retical curves are seen to increase rather sharply near Alach 
number 1.0. For buckled panels the increase is due to the 
use of steady -state linearized air forces which become infinite 
at M— 1. For flat panels the increase is thought to be 
associated with a change in flutter mode and decreased aero- 
dynamic damping. The curve for flat panels would have 
a finite ordinate at M— 1.0. 

As can be seen from figure 13, buckled panels have been 
fluttered up to a Afach number of 3. Flat panels were not 
fluttered over the same range because buckled panels ap- 
peared to be more susceptible to flutter, in general, and dur- 
ing a test it was difficult to prevent the thin, flat panels 
from buckling due to heating. (A temperature rise of 5° 
to 10° F was sufficient to induce buckling in many panels.) 

The agreement between theory and experiment for both 
flat and buckled panels is surprisingly good, inasmuch as the 
experiments were made on panels with a width-length ratio 
of 0.69, whereas the theories are for two-dimensional panels. 

CONCLUDING REMARKS 

A Rayleigh type analysis involving chosen modes of the 
panel as degrees of freedom has been used to treat the flutter 
of a two-dimensional flat panel supported at its leading and 
trailing edges and subjected to a middle-plane tensile force. 
The panel had a supersonic stream passing over its upper 
surface and still air below. The aerodynamic forces due to 
the supersonic stream were obtained from the theory for 
linearized two-dimensional unsteady flow and the forces due 
to the still air were obtained from acoustical theory. The 
still air beneath the panel was treated on the assumption 
that the still-air reservoir extended to infinity. Accordingly, 
once acoustic energy was radiated into this region, none of 
it was ever reflected. Such a situation is, of course, not the 
same as for a panel on a closed body but represents a first 


approximation for many practical cases. 

In order to study the effect of increasing the number of 
modes in the analysis, two and then four modes were em- 
ployed. The modes used were the first four natural modes 
of the panel in a vacuum with no tensile force acting. The 
analysis included the variables : Mach number M, structural 
damping, tensile force, density of the still air, and edge 
fixity (clamped and pinned). For certain combinations of 
these variables, stability boundaries were obtained which 
can be used to determine the panel thickness required to 
prevent flutter for any panel material and altitude. 

In contrast to some previous panel flutter investigations, 
the present results show that sufficiently thick panels are 
flutter free for the Mach numbers treated and suggest that 
this is true throughout the supersonic speed range. The 
low supersonic Mach numbers were found to be most critical 
from a design standpoint in the range examined (from 
M— 1.3 to M= 2.0). Tension was shown at M= 1.3 to have 
a marked favorable effect (also expected at all Mach num- 
bers) in reducing the thickness required to prevent flutter, 
and it was pointed out that one means of producing tension 
is by a static pressure difference between the upper and 
lower surfaces of the panel. Small amounts of structural 
damping were found to have a pronounced beneficial effect 
near and below M=-y[2 and essentially no effect at M= 1.56. 
In the neighborhood of M—^j 2 a small change in either 
ATach number or structural damping was found to cause an 
abrupt change in the thickness required to prevent flutter. 
At M—^2 a pinned-edge panel must be somewhat thicker 
than a clamped-edge panel in order to be flutter free: Still 
air below the panel was taken into account only at M= 1.3 
and was shown to have a moderate beneficial effect. For 
the still air is expected to have little effect because 
for this Mach number range the flutter mode is predomi- 
nantly the second natural mode, which radiates very little 
energy into the still air. 

The theories of the present report for flat panels and of 
Isaacs for buckled panels were compared with a few experi- 
mental results on panels clamped at leading and trailing 
edges over the Alach number range 1.2 to 3.0. The agree- 
ment was surprisingly good inasmuch as the experiments 
were made on panels with a width-length ratio of 0.69, 
while the theories are for two-dimensional panels. Over 
the Alach number range of the experiments it was found 
that buckled panels had to be thicker than flat panels in 
order not to flutter. The effect of restraining flat or budded 
panels on all four edges has not been investigated. Such 
restraint together with variation of width-length ratio will 
probably have a significant effect on the thickness required 
to prevent flutter. Another factor which requires investi- 
gation is built-in curvature of the panel in the streamwise 
or cross-stream direction. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., April 20, 1955. 
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APPENDIX A 

SOLUTION BY MEANS OF LAPLACE TRANSFORMS 

In a recent paper (ref. 6), a procedure is outlined for obtaining by means of Laplace transforms the exact solution 
for the flutter of a two-dimensional membrane which is subjected to a supersonic stream on one side and stagnant 
air on the other. This solution is called exact, inasmuch as the equation of motion for the system is solved directly 
without any limitation being imposed on the mode shape or frequency of flutter. Reference 6 also mentions that pure 
bending of a plate and the more general case in which plate bending and membrane stretching are combined could be treated 
in the same manner. The present report treats the latter case; namely, the flutter of a panel (plate) acted on by a middle- 
plane or axial force, such as tension, or compression less than the buckling load. In the body of the report this problem is 
solved by the generalized-coordinate approach, and the coordinates used are the normal modes of the panel with no axial force 
acting. In this appendix, for the sake of completeness, the solution to the same problem is derived by means of Laplace 
transforms to the point where numerical calculations can be made. The feasibility of applying the Laplace transform solution 
is examined, but no numerical results are obtained. 

The intergrodifferential equation to be solved is given by equation (2) which, upon substitution of the expression for 
p(x,t) obtained from equations (5) to (9), may be written as 

p Z”"-o?m A Z~ ^ [i^ 1 »,(«)/,(*- =0 (Al) 

or alternatively as 

~ t Z""-Jm A Z-f 2 Z"+^f^ [^(O^W+JJ /„(*-{) (^+»2fc) (A2) 

where 

w u (x)=Z'(x)-\- i 2 kZ(x) 

w,(x)=i2kZ(r) 

Iu( x )=e~ iz ’ z Jo (jjr xj 

I l (x)=H 0 ^(2kM\x \ ) 

Dividing equation (A2) by yields 


where 



The quantity u, in the formulas for a and 8 is the first natural frequency of the plate vibrating in a vacuum with no axial force 
F acting and Ki is the associated eigenvalue. (See table I.) In the case of the membrane, would be the first natural fre- 

/r 

quency of the membrane, a would be zero (D is negligible for the membrane), and / would be r 


Applying the Laplace transform 


L { Z{x ) } —Z((s ) — 


j; 


x Z{x)dx 
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to equation (A3) (using in the process pairs 4 and 7, p. 294, ref. 17) yields the transformed problem in the form 

a( s -'Z-^z Q -s 2 z l -sz 2 ~z a )~S(s 2 Z-.,z 0 -z l )-Z+~ Q [(s 2 +i4k.s-4k 2 )Z- (»+i2k)z 0 } I„(«)+g j 4 k 2 L ^ £ Z($)/,(x-$)$}-)=0 

(A4) 


where 2 0 =Z(0), z { — Z'(0), 2 2 =Z"( 0 ), z 3 =Z'"( 0), and J a (s) is the Laplace transform of I u (x). The Laplace transform in 
equation (A4) involving p„/p as a multiplier is the contribution of the perturbation pressure on the lower surface of the panel. 
Unfortunately, this transform does not appear to be obtainable here where the deflection Z is unknown. In the body of the 
report the effect of including the air below the panel is found to be moderate at a Mach number of 1.3 and reasons are given 
why this effect is expected to be even smaller at the higher Mach numbers investigated. In view of these facts and in view of 
the difficulty of handling the lower-surface term in equation (A4), this term will be omitted in the rest of this appendix — that 
is, treated as if p u were zero. 

Equation (A4) can therefore be reduced to 

[(os' 1 — 5s 2 — 1)4- e(s+'(2t) a 7„(6')]Z(s)=a(s-7,+6-.2 2 +2 3 )--5i', (A5) 

where e=l /4k 2 pf3 and has been dropped because it is zero for the present boundary-value problems. Thus the integro- 
differential equation (Al) has been reduced to the algebraic equation (A5). 

Now by means of pair 11, p. 294, and pair 55, p. 298, of reference 17 there is obtained 


Ai( s )— (s+L>)'+( 


Therefore, from equation (A5), after some algebraic manipulation 


where 


z <’>=W) + m m 

Q(s)=(os i —8s-—])' 2 + J— t-(s+i2ky 


[(*+L0M-(~) ] 


[«(,s ,2 2i 8 Si d - S;t) — 


-l/(.«0=(a.s 4 — 8ft - — 1) I 

X(s)= — f(s+i2k) 2 ^(s+w) 2 +(j0 J Hs^j+S^+Sa)— 8z t 
In polynomial form the quantities Q, M, and N are 


(Ah) 


(A7) 


Q(s)=^ l q r s U) - T 


AI(s)--- z,'y',m r a) s s r +2 3 Xl'Wr <3) s 6 r 

r=0 r=0 r=0 


(AS) 

(AO) 


xV(s)=SlXl«'r a) S G r + <:!) S 5_r + r 

r=0 r - 0 r=0 


(A 10) 



20 REPORT 1280 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 

The coefficients of the various series are given in the following table : 
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The exact inverse transform of equation (A7) requires the determination of the roots of Q(s ) (eq. (A8)). Since Q(s) is a 
tenth-order polynomial, its roots can be solved only approximately for specific values of the coefficients q T . An alternative 
procedure is to expand the quantity [(?(s)] -1 in a Maclaurin’s series (a procedure used in ref. 18), with the result that it may 
be expressed in the form 


L-=J 

Q(s) s w £i>s" 


(All) 


where 

<JoTq=1 


in 


<Z0 T n — — '^iQrTn-.r 
r = 1 


( n =i) 


and T with a negative subscript is to be interpreted as zero. 

When the series expansion for [<?(s)] _1 (eq. (All)) is substituted into equation (A7), the transform Z(s) becomes the sum 
of i nfin ite series with terms of the two distinct types 

A 

s m 

and 

« ft\ 


where m is a positive integer. The inverse Laplace transform of the first type of term is (sec pair 3, p. 295 of ref. 17) 


and of the second is (see pair 7, p. 294 of ref. 17) 



(»i— i)! 




{>-} 




(A12) 


(A13) 


where /„ (. r ) is defined following equation (A2). 
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Upon substituting equations (A9), (AlO), and (All) into equation (A7) and using equations (Al2) and (A13) to obtain 
the inverse transform of the resultant expression, Z(x) is given by 


where 


Z(x) = h i (x)z 1 +h 2 (x)s 2 +h 3 (x)s 3 (A14) 

« 1 T m (2 1 r n+r + 2 „ 5 rp (2) p x 

S +S S I (*-8" +r+4 J.ft)cfc 

«=o r=o «= u r=o (a+r-f-4)! Jo 

<n « r f m (3) r »+r+3 „ I 7' (3) rx 

*•«=»?. S fesr +S § sfess X 

In deriving equation (Al4) only one boundary condition— namely, 2 o ==Z(0) = 0 — has been used thus far. In order to 
obtain the solution for a plate restrained in a particular manner, it is necessary to impose three additional boundary conditions. 
These additional conditions for the plate with pinned and clamped edges are given in equations (3) and (4), respectively. 
By their use, one of the terms of equation (A14) is eliminated and two homogeneous equations in the two remaining unknown 
z t ’s are obtained. The borderline condition of harmonic oscillation, or the point at which flutter occurs, is obtained by setting 


5 equal to 

zero. 

Thus, the flutter determinant for the 

pinned-edge plate is 

hi (1) 

Ih( 1) 



K"(\) 


=0 

(A 15) 

V'( 1) 


j A-2 (1 ) 

IhO) 



k'(u 


=0 

(A16) 


and for the clamped-edge plate is 


where the determinant elements are given by 

£ (tSSoi+S, & ;>. J> 

- £ („V£ S (ZfrC >!+£, § J> -*>• 1 '«'•«>* 

4,(1 )- S § <rf^Ffjr + S § (« +;+4)! X 

s (SS ).+§ § x<->-™ 


A 3 (1)=S 2 


T„/h, ,3} 


» I 7’ ,, <3) 


,, (3) Cl 


(£)</£ 


«=<> r =0 (w+/*+ 8)1 w=0 r=o(ft + 

s (5S&+S § /> -« 

Each of the preceding elements contains integrals of the form 

/„ca/,s)= r 

which can be written in terms of the Schwarz functions /*(A/,«) (see ref. 9 or eq. (24)) as 


(A 17) 


(A 18) 
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where 

4= Jo ?I.(Qd£ 

Examination of the series in the elements ol' equations (A15) and (AlG) reveals that X of equation (A18) ranges at least 
between 0 and n-\-7 and at most between 0 and n + 9 . In order to obtain accuracy to four significant figures, at least the 
first eight terms of each series and the consequent ranging of X between 0 and 1G are probably required. Inasmuch as the 
Schwarz functions A have been tabulated for only the first few values of X, the use of equation (A18) would require the deter- 
mination of a rather extensive series of f\'x. An alternative and perhaps more efficient procedure would be to evaluate 
directly the integrals /,„ as given in equation (A 17) rather than to resort to the expanded form in equation (A18). 

Attention will now be given to the solution of the determinantal equations (Al. r >) and (AlG). A method of solution for 
parameters that were sought in the generalized-coordinate approach of the body of the report (that is, \/n and 
2k l =2k ( /{./’.it ) 1 «) will be outlined here. 

The elements of equations (AlG) and (A 1 (i ) are complex functions of the five parameters M, Ic, U (with £ = ()), J, and I//*. 
The most difficult parts of these elements to evaluate are the integrals generieally represented by /„, in equation (A 17), which 
are functions of the parameters M and k. Therefore, a convenient method of solution would be to fix the parameters M 
and k and preferably / and vary the remaining parameters 12 and 1/a in the left-hand side (hereinafter referred to as A) of 
equation (A15) or of equation (AlG). By varying If and 1/a over sufficiently broad ranges, an indefinitely large number of 
combinations of 12 and 1/a which cause A to vanish could be found. Each combination would define a point on separate sta- 
bility boundaries, such as those shown in figure 3. Each boundary could then be determined as completely as desired by 
varying k over a sufficient range and repeating for each chosen value of k the process of finding combinations of 12 and 1/a 
which cause A to vanish. 

As can be surmised, the numerical calculations would be extremely lengthy even apart from two other questions which 
arise; namely, which is the stable side of each boundary, and has the critical boundary been found which separates stable 
and unstable regions and thereby defines the thinnest panel that is stable? In the present report, therefore, the stability 
boundaries shown in figures 2 to 9 were calculated exclusively on the basis of the generalized-coordinate or modal approach. 
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APPENDIX B 


VELOCITY POTENTIAL AND RELATED INTEGRAL P n 


VELOCITY POTENTIAL 

The velocity potential 4n given in equation (9), which ap- 
plies to the lower surface of the two-dimensional panel 
shown in figure 1, will now be derived. The system consists 
of a panel of width c, which is part of an otherwise rigid 
surface of infinite ex ten l, oscillating harmonically with 
stationary air extending to infinity below. Thus, over 
the panel the normal velocity on the lower surface is 
Wi=iu>Z(x)e iul , while over the rest of the plane Wi—0. 

According to reference 10 the solution to this problem can 
be obtained from 


* = ~2 UJ W,, V ' IS 


(Bl) 


where w, is the given normal velocity at the element of area 
dS of the plane and <j> is the velocity potential at a point P 
which is at a distance r from <IS. From equation (Rl) the 
velocity potential at the surface of the panel may bo obtained, 
in terms of the coordinates of figure 1, as 


\ Z{m\ dy 

2 *■ do J— A U—ZY+y 2 ■' 


(B2) 


Upon making the substitution y=\x— £j cosh 6, the integral 
with respect to ?/ in equation (B2) may be written in the 
form 


* — V £) 2 +// 2 


/(*-«)= f“ — - ™ dy= f 

(.r-+) 5 +?r J 


-i cosh 1) 

= I e 0 de (B3) 


By means of equation (11) on page 180 of reference 19, 
equation (B3) becomes 

I {■ r ~£)= — >irHo r2) (B4) 

Substitution of equation (B4) into equation (B2) yields 

(B5) 

If the coordinates x and £ are nondimensionalized by dividing 
by the panel chord c, the form for <j>, given in equation (9) 
is obtained. 

INTEGRAL P„ 

The third term on the right-hand side of equation (22), 
namely, 

?„(*)=- 1 f ‘ ZJ&) log, {kM\x-l\)dt (Bfi) 

P TCJtt 

contains the singularity of the Ilankel function in equation 

( 20 ). 


As a first step in the evaluation of the improper integral 
in equation (B6), let 


and 


where 


1=9(1 — cos f) 


J-=2 (1 -cos '/') 


Z„(?)=Z 'jSm sin TOf 

m = 1 


(B7) 


(B8) 


2 

S m = Z„(i) sin mi <ii 

7T J (i 

In terms of equations (B7) and (BS), equation (BO) becomes 




Po J 
p 2 tt 




(B9) 


when 


<“ ('r /IrAf \ 

S,„J n sin i sin m{/oy r ^--^- |cos f —cos p\ jdi (BIO) 

Taking the derivative of equation (BIO) with respect to 
i P and making use of reference 20 to evaluate the resulting 
improper integrals yields 


-m = l Jo 


dp 


cos (m — 1) t — cos (m + 1) t . , , 

7 ) sin p di; 

cos f— cos p 


=oIjS» [sin(m — l)p— sin (m + 1)^] 

£ m = 1 


Integration of equation (Ii 11) gives 


(Bll) 


eos 2^+gS«[- 


os(m + l)^ cos(m — 


m + 1 


m- 


(Bl 2) 


The integration constant K in equation (B12) is determined 
by sotting p equal to x/2 in equations 1B10) and (B12) and 
equating the two resultant expressions. By so doing, it is 
found that 

K=lS x !og e ^ (B13) 

By means of equations (Bl2) and (Bl3), equation (B9) 
becomes 




cos 2 *)& + 
"cos (m-j-l)p cos (m — I 


(BI4) 

m =2 L m + 1 m — 1 J J v 
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The mode shapes Z n in equation (B6) are now approxi- 
mated by the finite sine series 

10 

Z n {£) ~^jS T sin rx (B 1 5) 

T—l 


The constants S m in equation (B14) are obtained from the 
expression following equation (B8), with the result that 
equation (B14) can be written as 


£(*)=- I (*K 

P 4 r= i 

where 

r,(j')=log e cos 


(BIG) 


cos (r+1 )i- - cos (r- 1 )+ (r ^ 2) 


The form given in equation (BIG) was used to obtain equa- 
tion (28). Of interest is the fact that only Li depends on Jc 
and M. The term P„ is therefore comparatively simple to 
include in equation (22). 

The coefficients s r in equations (Bl5) and (B16) for the first 
four modes of the plate with clamped edges are given in the 
following table: 



Mode 1 

Mode 2 

Mode 2 

Mode 4 


0. 06T>13 

0 

-0.22917 

0 

* 2 

0 

0. 04119 

0 

0. 40890 

•'Ll 

-0. 29502 

0 

-0. 50807 

0 


0 

-0. 40801 

0 

0. 27852 

nr, 

0. 029704 

0 

0. 01 124 

0 

«rj 

0 

0.098515 

0 

-0. 01524 

H ; 

0. 001215 

0 

-0. 18740 

0 

.Ss 

0 

-0.002411 

0 

(1. 2550(1 

HU 

0. 000290 

0 

0.021097 

0 

Kid 

0 

0.002858 

0 

-0. 028044 


and, as in equation (B7), 

x=l (1 —cos i) 


The coefficients .v 1( > for the second and fourth modes were 
obtained bv forcing the slope of Z n , as given by equation 
(Blf>), to be zero at ^=0 (that is, at x=0). A similar table 
can be easily calculated for the pinned-edge plate. 
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